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Abstract: The ephemeris data format of legacy GPS receivers is improper for positioning 
stationary pseudolites on the ground. Therefore, to utilize pseudolites for navigation, GPS 
receivers must be modified so that they can handle the modified data formats of the 
pseudolites. Because of this problem, the practical use of pseudolites has so far been 
limited. This paper proposes a pseudolite -based positioning system that can be used with 
unmodified legacy GPS receivers. In the proposed system, pseudolites transmit simulated 
GPS signals. The signals use standard GPS ephemeris data format and contain ephemeris 
data of simulated GPS satellites, not those of pseudolites. The use of the standard format 
enables the GPS receiver to process pseudolite signals without any modification. However, 
the position output of the GPS receiver is not the correct position in this system, because 
there are additional signal delays from each pseudolite to the receiver. A post-calculation 
process was added to obtain the correct receiver position using GPS receiver output. This 
re-estimation is possible because it is based on known information about the simulated 
signals, pseudolites, and positioning process of the GPS receiver. Simulations using 
generated data and live GPS data are conducted for various geometries to verify the 
proposed system. The test results show that the proposed system provides the desired user 
position using pseudolite signals without requiring any modifications to the legacy GPS 
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receiver. In this initial study, a pseudolite-only indoor system was assumed. However, it 
can be expanded to a GPS-pseudolite system outdoors. 

Keywords: GPS; pseudolite; indoor navigation 



1. Introduction 

Over the past few decades, there have been a number of studies on Global Navigation Satellite 
Systems (GNSS). Many people around the world nowadays utilize GNSS receivers on mobile devices 
for various kinds of applications. With the spread of location-based services and applications, the need 
for position information has grown even in places where GNSS signals are not available, such as 
indoors, in mountainous or urban canyons, and in other GNSS shadow areas. Many studies have been 
done to increase the availability of GNSS in these challenging environments [1-3]. A pseudolite is a 
GNSS -like signal generator and transmitter that can be used as an additional signal source to augment 
GNSS. Using pseudolites has many advantages over other methods [4-6]. Employing multiple 
pseudolites, it is even possible to develop an independent positioning system where the GNSS signals 
do not reach. However, there are some drawbacks to using pseudolites in positioning systems, 
including the following [7]: 

(1) Near/far problem 

(2) Phase noise 

(3) Receiver modifications 

This study focuses on the receiver modification problem. The receiver modification problem 
refers to the need to modify GNSS receiver hardware or firmware to acquire accurate data using a 
pseudolite-based positioning system. In a pseudolite-based system, only GNSS receivers developed for 
special uses can calculate their own positions, because the pseudolite positions cannot be delivered in 
ephemeris format. Therefore, legacy GNSS receivers such as those in smartphones or mobile devices 
cannot be used with pseudolite signals. 

To overcome this problem, a study was done that suggested a modified navigation message to 
deliver a pseudolite's position information [8]. The modification is simple yet tricky. Pseudolites 
'cheat' the GPS receivers by filling in some ephemeris parameters with particular numbers to make the 
resulting coordinates independent from time and fixed at required values. However, the scale factor 
and the number of data bits used in GPS ephemeris data format is limited [9]. Therefore, it is 
impossible to make the pseudolite position perfectly independent from time, and the accuracy and 
precision of the pseudolite position remains limited. 

This paper proposes a new positioning system for pseudolite signals that does not require any 
modification of legacy GPS receivers. The proposed system uses pseudolites which generate each 
channel signal of the GPS simulator and legacy GPS receiver without any modification. The user 
position is estimated using the output of the legacy GPS receiver and the information about pseudolites 
and simulated GPS satellites. Section 2 shows the limits of expressing pseudolite positions in GPS 
ephemeris format. Section 3 proposes a new positioning algorithm using pseudolite signals. The 
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system concept and detailed algorithms are described. Simulation results will be presented in Section 4. 
The test was performed using simulated measurements and live GPS LI data with additional manual 
delays. The results show that an accurate position can be obtained using the proposed algorithm 
without any receiver modification. For convenience, the GPS LI signal with C/A code was assumed 
during the study. 

2. Ephemeris Parameter for Pseudolites 

Before switching to a pseudolite -based positioning system, it makes sense to check for the 
possibility of expressing pseudolite positions using GPS ephemeris parameters. Table 1 gives the 
definitions of the orbital parameters using typical Keplerian orbital parameter terminology, and it also 
lists the relevant ephemeris parameter types. Legacy GPS receivers determine GPS satellites' orbits 
using these parameters [9]. The parameters are two's complements with the sign bit (+ or -) occupying 
the MSB. Exceptions to this rule are e , A^''^ , and t^^ . 



Table 1. Selected ephemeris parameter of GPS LI signal. 



Ephemeris 
Data 


Definition 


No. of Bits 


Scale Factor 
(LSB) 


Units 


M, 


Mean anomaly at reference time 


32 


T'' 


Semi-circles 


An 


Mean motion difference from computed value 


16 




Semi-circles/sec 


e 


Eccentricity 


32 




Dimensionless 


^1/2 


Square root of the semi-major axis 


32 




Meters 


Q, 


Longitude of ascending node of orbit plane at weekly epoch 


32 




Semi-circles 




Inclination angle at reference time 


32 


2-'' 


Semi-circles 


CD 


Argument of perigee 


32 


2-'' 


Semi-circles 


Q 


Rate of right ascension 


24 


2-'' 


Semi-circles/sec 


/ 


Rate of inclination angle 


14 


2-'' 


Semi-circles/sec 



The number of bits and scale factors of GPS ephemeris parameter are designed for GPS 
satellites with radii of approximately 26,000 km. In order to employ GPS ephemeris format for 
pseudolites located on Earth, Q should be identical to the rotation rate of the Earth, which is 
7.2921 151467e-5 rad/s [9]. Because this is out of the range defined for Q , which has 24 bits and a 
scale factor of 2x10-43, it is impossible to express accurate pseudolite positions on the Earth using 
GPS ephemeris data format. Therefore, the ephemeris format for GPS is not suitable for pseudolites on 
the ground. However, legacy GPS receiver can only process the ephemeris format of GPS. Therefore a 
new method is needed to utilize pseudolites with legacy GPS receivers. 

3. Pseudolite-Based Positioning System for Legacy GNSS Receivers 

In this section, the proposed pseudolite-based positioning system will be described. 

3.1. System Configuration 

The proposed system is composed of four main parts: the pseudolites, the data server, the reference 
station, and the user terminal. Figure 1 illustrates the total system configuration. 
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Figure 1. Illustration of proposed system. 
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3.1.1. Pseudolites 

In this system, pseudolites generate simulated GPS LI signals and transmit them for users. The 
necessary information for the simulation is provided by the Data Server. Using the simulation scenario, 
pseudolites generate one channel signal of the simulation scenario. Figure 2 shows the structure of the 
pseudolites used in this system. The CPU controls overall components of the pseudolite according to 
the simulation scenario received from the data server. The Digitally Controlled Oscillators (DCO) is 
the clock source of the pseudolite. By providing control inputs to the DCO, the timing of the 
transmission of the pseudolite signal can be controlled. Also, the timing of the signals of the plural 
pseudolites can be synchronized by providing proper control inputs to each pseudolite. Generation of 
the clock control input for the synchronization of the pseudolites will be described later. 

Figure 2. Pseudolite structure. 
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As described, the pseudolites in this system do not provide any information about themselves 
(pseudohtes) in their signals. They just broadcast simulated GPS signals as there were no pseudolites. 
The reason of using this kind of pseudolites is to make the legacy GPS receiver process the pseudolite 
signals without any modifications. 

3.1.2. Reference Station 

The reference station whose position is accurately known is needed for the clock synchronization of 
the pseudolites. It includes a reference GPS receiver to collect measurements of the pseudolite signals. 
The reference station provides the data server with these measurements. Then the data server generates 
the clock control inputs for each pseudolite. 

3.1.3. Data Server 



The data server plays two main roles in this system. One is to provide the pseudolites and the 
user terminal with necessary information. It gives the information about the simulation scenario to the 
pseudolites so that each pseudolite can generate a simulated GPS signal. Also, it provides the 
user terminal with the information about simulated GPS satellites and pseudolites. The delivered 
information include the positions of pseudolites and simulated GPS satellites which are required to the 
post-calculation module. It also includes the information about the assignment of the pseudolites and 
corresponding GPS satellites. Using this information, the user terminal distinguish the pseudolites 
using PRNs of corresponding GPS signals. The data link between the data server and the user terminal 
can be WiFi or any other possible means. 

Figure 3. Control Block Diagram of Pseudolite Clock Synchronization [10]. 
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The other role of the data server is generating clock control for the synchronization of the 
pseudolites. However, this paper will not describe the synchronization of the pseudolites in detatil, 
because the synchronization can be achieved by using an existing method. This paper follows the 
method using a reference station receiver to synchronize the clocks of pseudolites [10]. Figure 3 shows 
the control block diagram of the pseudolite clock synchronization. According to the paper, the relative 
clock error of the pseudolites can be calculated by differencing the carrier phase measurements of the 
master pseudolite and a slave pseudolite, and subtracting the known distance values. Then the properly 
designed loop filter F(s) estimates the control value of each pseudolite clock. The output of the DDS 
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(Direct Digital Synthesizing) is then controlled to eliminate the clock error a of the TCXO of the slave 
pseudolite. The experiment results showed cm-level accuracy in CDGPS positioning based on the 
clock synchronized pseudolites. 

3.1.4. User Terminal 

The user terminal includes a general legacy GPS receiver, a data communication module, and a 
post-calculation module. The GPS receiver calculates a position using pseudolite signals. However, 
this a priori position is not correct, because the simulated GPS signal includes additional delays due 
to the use of pseudolites while the receiver does not know the existence of the pseudolites. The 
post-calculation module re-processes this false output of the GPS receiver using information about 
GPS satellite positions and pseudolite positions. The detailed algorithm of the post-calculation module 
will be described later. Figure 4 illustrates the structure of the user terminal. A smartphone is a good 
example of a user terminal, as it includes a GPS receiver and data link. A platform to adapt the 
post-calculation algorithm is also prepared. 

Figure 4. Block diagram of user terminal. 
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3.2. User Algorithm Concept 

In this system, pseudolites simulate GPS signals received at a particular receiving point , and 
each pseudolite corresponds to the simulated GPS satellite. If the distances between the user terminal 
and each pseudolite are identical, as in Figure 5a, the output of the GPS receiver in the user terminal 
becomes the receiving point . The GPS receiver regards the received signals transmitted from not 
the pseudolites but the GPS satellites, as shown in Figure 5b, and the additional signal delay due to 
the use of pseudolites becomes an additional common bias that does not affect the position result. 

If the distances between the user terminal and each pseudolite are different from each other, as in 
Figure 6a, the output of the GPS receiver in the user terminal will also differ from the receiving point. 
From the perspective of the GPS receiver, the signals from Satellite 1 are delayed by d^+d^, and 
signals from Satellite 2 are delayed by d^+d^. Therefore, the output of the GPS receiver is neither the 
receiving point R^ nor the user position. 
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Although the output of the GPS receiver is not the desired one, it contains some information about 
the delays , and J2 • Therefore, using the output position of the GPS receiver and the information 
about simulated GPS satelhtes and pseudolites obtained from the data server, it is possible to estimate 

and ^2 ^iid use this information to re-estimate the desired user position. 



Figure 5. (a) A pseudolite-based navigation system with identical travel time for both 
signals (2-D case); (b) Position result as given by the GPS receiver in the user terminal 
(2-D case). 
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Figure 6. (a) A pseudolite-based navigation system with different travel times for each 
signal (2-D case); (b) Position result given by a GPS receiver in the user terminal (2-D case). 
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3,3, Pseudolite-Based Positioning Algorithm for Legacy GNSS Receivers 

In this section, the proposed user algorithm will be described with detailed formulas. Figure 7 
illustrates the j-th pseudolite and user terminal. 
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Figure 7. Illustration of j-th pseudolite and user terminal (Simulated and Real Simuation). 
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R- : position vector of the j-th pseudolite 

: position vector of the user 
d^- : distance between the user and the j - th pseudolite 
e^- : unit vector from the user to the j - th pseudolite 
R^ : position vector of the GPS satellite simulated by the j-th pseudolite 
R^ : position vector of the receiving point for the simulated signal 

: distance between the receiving point and the GPS satellite simulated by the j-th pseudolite 
ej. : unit vector from the receiving point to the GPS satellite simulated by the j-th pseudolite 

3.3.1. Positioning Process of GPS Receiver 

The legacy GPS receivers in the user terminal estimates the user position using pseudolite signals. 
The pseudorange measurement for the j-th pseudolite can be expressed as Equation (1): 



(1) 



where T is the simulated tropospheric delay of the j - th GPS signal, is the simulated ionospheric 
delay of the j-th GPS signal, and is the tropospheric delay of the j-th pseudolite signal. B is 
the clock bias of the GPS receiver, and s^^ is the thermal noise. This study assumes that , , and 

are either negligible or are eliminated by the model. The pseudorange measurement model is thus 
simplified into Equation (2): 



(2) 



Equation (2) is the measurement model for true situation. However, the GPS receiver does not 
know the existence of pseudolite. Therefore the measurement model used by the GPS receiver will 
follow the general pseudorange measurement model as Equation (3): 



r^rcv u,rcv rev nJ 



R'-R 



(3) 
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d^^^^ is the distance between the j-th GPS sateUite and the user, B^^^ is the clock bias of the 
receiver, and s , is thermal noise. If a priori estimations of the user position and clock bias are 

Prcv 

established so that x^^ = B^^ j , the estimated pseudorange measurement should be as follows: 



Prcv 



Hrcv 



(4) 



The pseudorange measurement residual can be calculated by subtracting Equation (4) from 
Equation (3): 



Spi =ni -pi =Rj-R -Rj-R 

r^rcv rrcv r rev u,rcv u,r 

Linearization is performed at current estimation: 



+ B_-B_+s , -s.. 



rev rev 



(5) 



(K.rcv-K,rcv) + 



,rcv ,rc 



(5.cv-5.cv) + (V ) 



(6) 



For m satellites, m equations are determined in the same way, and a matrix equation is constructed: 



5P, 



-e' 1 

u,rcv 



-e"' 1 

u,rcv 



5e 



Prcv 



5s 



Prcv 



(7) 



5p =H dx +ds 

r^rcv rev rev p^^ 

The error state can be estimated by the Least Square Solution [11]: 

5x ^\5R ^ 5b\ ^{h Vh ^5p 

rev u,rev \ rev rev J rev r^rc 

Next, state variables are updated: 



X = X +5x 

rev rev, old rev 
^ rev, old "^rev 



(8) 



(9) 



These processes are repeated until the state converges. The converged state may be expressed as 



Kev,evg = \_K,rev,evg^ ^rev,evg ] * ^hls Is thc output of thc GPS recclvcr aud the a priori position for the 
post-calculation module. 

3.3.2. Estimation of Measurements from GPS Receiver Output 

The post-calculation module have to estimate the user position using the output of the legacy GPS 
receiver. If the GPS receiver in the user terminal provides pseudorange measurements, the receiver 
output position can be ignored and the desired user position can be calculated using the output 
measurement. However, some legacy GPS receivers and some mobile devices that are equipped with 
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GPS receivers does not provide these measurements to the user. In these cases, the original 
measurements must be estimated from the output position of the GPS receiver. 

As described in Section 3.2, though the GPS receiver output calculated in Section 3.3.1 does not 
give the correct position, it contains some information about the original measurements. The 
measurements from the position output can be estimated via Equation (10). Satellite position W is 
provided by the data server: 



R'-R 



u,rcv,cvg 



+ B 



rcv,cvg 



(10) 



3.3.3. Final Estimation of the User Position 

Using given pseudorange measurements, finally the post-calculation module can estimate the 
user position. From Equation (2), pseudorange measurements can be rewritten using position vectors 
as in the following equation: 



pJ =\RJ-R\ + 



R^-R..\ + B + e 



(11) 



Again, an a priori estimation of the state vector as jc = |^j^y^ 5 J is assumed. The estimated 
pseudorange measurements may then be described by Equation (12): 



p^=\r^-r\+ 



R^-R 



+ B + £^ 



(12) 



By subtracting Equation (12) from Equation (11), the pseudorange measurement residual is 
calculated using the following equation: 



Sp' = p' - p^ 



R'-R 



R'-R 



+ B-B + S . 

p' p' 



(13) 



Equation (13) is linearized at the a priori estimation: 



Sp 



J ^ dSp' 



d5p' 



dB 



(B-B) + 



I SB 



(14) 



P' 



The same equations can be used for m pseudolites: 



'5p^~ 




'-el I 




Ss 1 

p 
















I > 












dp"" 




.-^f 1. 








Sx 




dp 


V 

H 




5sp 



(15) 



5p = HSx + SSp 
The error state is estimated by Least Square Solution [11]: 

5x = [SRj SbJ = {h^h)~' H^Sp 



(16) 
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The original state is updated using estimated error state: 



(17) 



The final estimation is obtained by repeating these processes until the states converge. The flow 
chart of the algorithm for the post-calculation module of the user terminal is illustrated in Figure 8. 

Figure 8. Algorithm for the Final Estimation of the User Position. 
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In Figure 8, blocks in the dashed box is exactly the same as the position estimation process of 
general GPS receivers except the nonlinear equation used in Equation (12) (red box). The differences 
are coming from the existence of the pseudolites between the imaginary receiving point and the user 
receiver antenna. 
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4. Simulation Results and Discussion 

Unfortunately, it was hard to be equipped with pseudohtes which operates as a GPS simulator for 
this system. Also we could not make a hardware GPS simulator by ourselves. Therefore, the only way 
to verify the system was the test using simulated data. By the results, we tried to show the feasibility of 
the proposed new concept of the pseudolite based navigation system. 

Due to the reasons above, the simulated data and SDR (Software Defined Radio) are used to verify 
the proposed system. The SDR used in the verification is a platform to test various receiver algorithms. 
The operation of legacy GPS receivers also can be implemented in the SDR. In this paper, the SDR 
which is designed to operate just as the legacy GPS receivers is used for the test. Therefore, the use of 
the SDR instead of the legacy GPS receiver seems reasonable. 

The rest of this section shows some simulation results to verify the proposed system. First, a 
simulation is performed to show the feasibility of the algorithm. Second, the same test is carried out 
under poor geometries. Third, a test using live GPS signal is performed. Some manual delays due to 
the pseudolites are added to the live GPS measurements in this test. 

4.1. Feasibility Test 

The feasibility of the proposed algorithm is tested using five simulated satellites and five 
corresponding pseudolites. Figure 9 shows the procedure of the feasibility test. 

Figure 9. Simulation procedure. 
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First of all, true values are provided for the initial setting module. The positions of simulated GPS 
satellites and pseudolites are shown in Figure 10. For convenience, the receiving point is set at (0,0,0), 
and satellite positions are set as they were on the sphere of radius of 22,000 km centered at the 
receiving point. Local user positions were set at (0,0,0) and (1,1,0). PDOPs for GPS satellites and 
pseudolites were then calculated to be 1.7489 and 2.0468, respectively. 
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The initial setting module stores each true value to the proper variables. Using these values, 
pseudorange measurements are generated according to Equation (2). Tropospheric and ionospheric 
delays are assumed to be eliminated by the model or by other methods. 

Figure 10. (a) Skyplot of simulated GPS satellites viewed from the receiving point 
(b) Local positions of pseudolites and user. 
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Additional Gaussian random noise is added to the measurements. Variance of the measurement 
noise is set to 1 meter for all pseudolites. Figure 11 shows the generated pseudorange measurements 
for each pseudolite received by the User #2. 

Figure 11. Generated Pseudorange Measurements. 
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The GPS receiver calculates the a priori position using generated measurements as if it were the 
general legacy GPS receiver. The post-calculation module uses this a priori position, the true GPS 
satellite position, and the true pseudolite position to re-estimate the final position. Instead of the data 
server, necessary information are provided to the post-calculation module directly in this 
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simulation. Figure 12a shows that the a priori position of User #1 calculated by the algorithm stated in 
Section 3.3.1 is not biased. In this case, the distances from all pseudolites to the user are the same. 
Therefore, all pseudolite signals are subject to the same delay, which affects only the clock error, not 
the position measured. The horizontal and vertical drms errors are calculated to be 1.2977 m and 
2.2927 m, respectively. Unlike User #1, User #2 got biased a priori positions from the GPS receiver, 
as shown in Figure 12b. The distances from the pseudolites to the user are different from each other in 
this case, and they affect the position results. However, using the proposed algorithm, the final position 
is estimated with reasonable accuracy. The horizontal and vertical drms errors are calculated to be 
1.2977 m and 2.2927 m, respectively. The a priori position of User #2 has much bias error than that of 
User #1. However the error variance is calculated as the same value. 

Figure 12. (a) a priori position for user #1 (b) a priori position for user #2. 
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As described in Section 3.3.2, though the GPS receiver output does not give the correct position, it 
contains some information about the original measurements. The measurements from the position 
output can be estimated via Equation (10). Figure 13 shows the difference between the originally 
generated and re-estimated pseudorange measurements. The residual values are affected by the 
geometry of the simulated GPS satellites and pseudolites. 

Figure 13. Estimated Pseudorange Measurements. 



Pseudorange estimation error 




Sensors 2014, 14 



6118 



Using the a priori positions and the estimated pseudorange measurements, the final estimation of 
user position can be calculated as described in Section 3.3.3. Figure 14 shows the a priori and final 
position of User #1 and User #2. In the final estimation, the bias error is eliminated for not only the 
User #1, but also for User #2. The values of the horizontal and vertical drms and bias errors became 
smaller than those of the a priori estimates. 

The results show that using the proposed system, user position can be estimated despite reliance on 
pseudolites combined with an unmodified legacy GPS receiver. The algorithm works regardless of the 
existence of bias in the output of the GPS receiver. However, it can be affected by the geometry 
of simulated GPS satellites and pseudolites. In the next section, simulation results for poor DOPs 
are presented. 
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Figure 14. (a) Estimation results for user #1 (b) Estimation results for user #2. 
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4,2, Poor Geometry Test 



Simulated GPS satellites never have severely poor geometry, because simulation allows the satellite 
configuration to be manipulated as necessary. However, the geometry may become less than ideal if 
the receiver unexpectedly loses a pseudolite signal corresponding to a GPS satellite. Therefore, the 
proposed algorithm is tested for cases involving poor GPS and pseudolite constellation geometry. 
Once again, five GPS satellites using five corresponding pseudolites are simulated. The simulation 
procedure and settings except the GPS satellites, pseudolite, and user positions are the same as the 
feasibility test. 

Figure 15 illustrates the results for good and bad GPS satellite geometries. For each case, the DOPs 
for GPS are 1.7489 and 14.5487, respectively. The pseudolite DOP is calculated to be 2.0468. With 
higher DOPs, the a priori position gets worse, and the final position error is also affected. However, 
unlike the biased a priori position, the final estimated position still matches the true position. 
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Figure 15. (a) Results for good DOP of GPS satellites (b) Results for poor DOP of GPS satellites. 
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To test for even worse cases, poor pseudolite geometries are used in addition to poor GPS 
geometry. Figure 16 illustrates the results for high pseudolite DOPs. In each case, the DOPs for the 
pseudolites are 14.1849 and 25.3421, respectively. For the same a priori positions, poor pseudolite 
geometry makes the final estimation even worse. Again, unlike the biased a priori positions, the final 
results match the true user positions. 

Figure 16. (a) Results for poor DOP of pseudolites (b) Results for worse DOP of pseudolites. 
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It can be assumed that any time a pseudolite system is developed, the pseudohtes will be arranged 
in an effective configuration. However, if the user's GPS receiver loses some pseudolite signals, the 
system will be left with poor geometry. Simulation results show that even when the geometry gets 
worse, the proposed algorithm works with only a slight loss of accuracy. 

4.3. Test Using Live GPS Signals 

The proposed algorithm is then tested using live GPS signals. Due to a lack of available 
GPS-simulating pseudolites, a live GPS LI signal is used instead of GPS -simulating pseudolites. Live 
GPS signals were collected using a GP 2015 RF board and an NI DAQ card [12]. The raw data was 
stored as a binary data file on a PC. The stored data is then processed by a Software Defined Radio 
(SDR) developed by the Seoul National University GNSS Laboratory [13]. The SDR output 
pseudorange includes the distance from the GPS satellites to the GPS antenna, clock errors, tropospheric 
and ionospheric delays, etc. The signal delay from the pseudolite to the user was manually generated 
and added. This simulated pseudorange was then used to test the proposed algorithm. Figure 17 
illustrates the simulation procedure. 
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Figure 17. Simulation procedure. 
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Four GPS satellites and four corresponding pseudolites are used in this test. The GPS and 
pseudolite constellation is shown in Figure 18. To test for different geometries, user positions are set to 
both (3,3,0) and (20,20,0) meters. The PDOP for pseudolites with respect to User #1 and User #2 are 
2.2371 and 5.7674, respectively. The PDOP for the GPS satellites is 15.0590. 

Test results for each case are illustrated in Figure 19. For both cases, the GPS receiver output is 
biased. The amount of bias is greater when the geometry is poor. However, the final estimated results 
place the users in the correct positions. The horizontal and vertical drms errors are 1.3709 m and 
1.7564 m for User #1. For User #2, horizontal and vertical drms errors are calculated to be 2.7491 m 
and 2.4726 m. 
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Figure 18. (a) Skyplot of live GPS satellites (b) Local positions of pseudolites and user. 
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Figure 19. (a) Estimation results for User #1 (good DOP) (b) Estimation results for User 
#2 (poor DOP). 
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5. Conclusions 



One of the major problems inherent in pseudolite-based positioning systems is the receiver 
modification problem. Due to this problem, pseudolites cannot be used with general legacy GPS 
receivers. To resolve this problem, this study proposed a pseudolite-based positioning system that can 
be used by unmodified legacy GPS receivers. The proposed system uses GPS -simulating pseudolites 
whose signal can be processed by legacy GPS receivers. The position output of the GPS receiver is 
then reprocessed using the position information of the simulated GPS satellites and pseudolites. 
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Several simulation tests were conducted using MATLAB to verify the proposed system. It was 
shown that the proposed system results in an accurate user position, unlike the biased GPS receiver 
outputs. Even under poor geometries of pseudolites or simulated GPS satellites, the proposed 
algorithm successfully obtained users' positions. Because GPS -simulating pseudolites were not 
available, live GPS signals were used to verify the system. Intentional delays were added to the 
pseudorange measurements in order to simulate the GPS signals transmitted by pseudolites and 
received by a local user's receiver. The results show that the algorithm works well for various 
pseudolite geometries and corresponding simulated GPS satellite geometries. 

This paper proposed a new concept of the pseudolite based navigation to overcome the receiver 
modification problem of the pseudolite-based positioning systems. The feasibility of the system is 
shown by test results using simulated data. In the future, proposed system should be tested by 
physically real data using real GPS-simulating pseduolites and general legacy GPS receivers. 

If this system concept matures, then pseudolite-based indoor or outdoor navigation will be available 
for many people with GPS receivers in their smartphones. The only other requirement would be an 
application that reprocesses the output of the GPS receiver. 
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